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ABSTRACT 

Observations from multiple gamma-ray telescopes have uncovered a high energy gamma-ray source 
spatially coincident with the Galactic center. Recently, a compelling model for the broad-band gamma- 
ray emission has been formulated which posits that high energy protons emanating from Sgr A* 
could produce gamma-rays through ir° decays resulting from inelastic collisions with the traversed 
interstellar gas in the region. Models of the gas distribution in the Galactic center region imply that 
the resulting 7-ray morphology would be observed as a point source with all current telescopes, but 
that the upcoming Cherenkov Telescope Array (CTA) may be able to detect an extended emission 
profile with an unmistakable morphology. Here, we critically evaluate this claim, employing a three 
dimensional gas distribution model and a detailed Monte Carlo simulation, and using the anticipated 
effective area and angular resolution of CTA. We find that the impressive angular resolution of CTA 
will be key to test hadronic emission models conclusively against, for example, point source or dark 
matter annihilation scenarios. We comment on the relevance of this result for searches for dark matter 
annihilation in the Galactic center region. 

Subject headings: (ISM:) cosmic rays — gamma rays: theory — gamma rays: observations 



1. INTRODUCTION 

Early observations from the High Energy Spectroscopic 
System (H.E.S.S) opened a new window into 7-ray obser- 
vations of the Galactic Center (GC) Q region, including 
the detection of a bright TeV point source localized to 
within 1' of the GC. The spectrum of the 7-ray source 
is fairly hard, following a power-law a = -2.2 +/- 0.09 
(stat) + / - 0.15 (sys) with a high -energy cutoff exceeding 
10 TeV (jAharonian et all |20M ). Further observations 
succeeded in localizing the center of the point-source to 
within 13" of the GC, and found 85% of the total 7- 
ray emission to be confined within 3 pc (1.2') of the GC 
(jAcero et al.ll2010D . 

While the observed morphological details strongly sug- 
gest that the TeV 7-ray signal stems from a point-source 
spatially coincident with the black hole at the GC, the 
steady-state nature of the H.E.S.S. source indicates the 
emission may be originating farther from the GC. While 
lower energy X-ray and radio observations have uncov- 
ered significant variability from Sgr A* on timescales 
stretching from minutes to years, no v ariability has yet 
been observed in 7-ray observations (jAharonian et al.l 
2009). Most notably, a simultaneous observation with 
H.E.S.S. and Chandra found that an X-ray outburst ob- 
served in 2007 wa s not correlated with any change in 
the 7-ray emission (jAharonian et~aL|[2008l ). This implies 
that the source of the 7-ray emission may be distinct 
from the source of low-energy photons. Several mod- 
els have been posited which would naturally explain an 
intense TeV 7-ray emission which is uncorrelated with 
the lower-energy regime, including photons from dark 

1 Throughout this work, we will employ the term Galactic Center 
(GC) to refer to both the dynamical center of the Milky Way, as 
well as to the position of the radio source Sgr A*, which we will 
consider to be equivalent. 



matter annihilat i on (jHooper et al.|[200l: iProfumol 120051: 
lAharonian et al.l 1200 6a). as well as 7r°-decay resulting 
from the emission of high energy protons from the central 
black h ole and their subsequent co l lisions with Gala c- 
tic gas (lAharonian fc Neronovl 120051: ILiu et al.1 12006allbl: 
iFrver et al.ll2007t iBallantvne et al.ll2007l ). 

With the launch of the Fermi-LAT in 2008, the win- 
dow was opened to observe the GeV 7-ray spectrum 
with simi lar angular and en ergy resolution to that of 
H.E.S.S. (jAtwood et al.ll2009D . unveiling a distinctly dif- 
ferent spectral shape from the very high-energy regime. 
Specifically, an excess of 1-10 GeV 7-rays was uncov- 
ered in the GC compared to that expected from an 
extrap olation of the TeV source spectrum to GeV en- 
ergies dVitale et al l 120091: iHooper fc Goodenoughl 120111: 
IHooper fc Linden! 120111: Hooper et al.l 120121) ! Several 
models have been postulated to explain this excess 
emission, including the a nnihilation of light, lept o philic 
dark matter particles ([Hooper fc Goodenoughl 120111 : 
IHooper fc Lindenl 120111) and emissio n from millisecond 
pulsar s (|Abazaiianll201lD . Recently, iChernvakova et all 
(|2011[) re-examined an extension of the hadronic scenario 
described above down to GeV energies and found that the 
entirety of the GeV- TeV spectrum could be explained by 
inelastic processes initiated by protons whose spectrum 
would follow a single power-law. The softening of the 7- 
ray signal at energies ^10-100 GeV was then enforced by 
fine-tuning the diffusion constant in order to produce dif- 
fusive propagation at GeV energies and rectilinear prop- 
agation at TeV energies. 

Subsequently, iLinden et al.l (|2012f ) examined the ex- 
pected mor phology of the had r onic e mission model de- 
scribed by IChernvakova et all ([201 ID . using a realistic 
model for the morphology of Galactic gas in the in- 
ner 1 pc around Sgr A* as determined by iFerrierel 
(|2012t ) . They found that the morphology of TeV emission 
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closely matched observations by Aharonian eFaLl (|2006al) 
signaling that 85% of the 7-ray emission from the GC 
was confined to wit hin ^3 pc of Sgr A*. Additionally, 
iLinden et all (|2012| ) found that the energy dependence 
of this morphology is minimal, and thus the majority of 
GeV emission detectable by the Fermi-LAT should also 
reside within 3 pc of the central black hole, which may 
be in tension with an o bserved extension of the GeV 
7-ray source at the GC (|Hooper fc Goodenoughl [20111 : 
iHooper fc Lindenl 120 111 ). We note that disentangling 
an extended emission in this region is especially prob- 
lematic, give n the poor knowled ge of the di ffuse Galac- 
tic emission (|Vitale et a l. 2009). Finally, ILinden et all 
(2012) noted that while the hadronic emission model 
should appear point-like to all current 7-ray instruments, 
the upcoming Cherenkov Telescope Array (CTA), should 
observe an extended spatial morphology which would dis- 
tinguish this scenario from other point-source emission 
mechanism occurring at the position of Sgr A*. 

In this letter, we closely investigate several models for 
the TeV 7-ray emission at the GC which may be ob- 
served by CTA, including a 7-ray point source at the po- 
sition of Sgr A*, p-p collisions due to hadronic emission 
from the pos i tion of Sgr A* (jChernyakoya et all 120111: 
Linden et al.l I2012D and dark matter (IHooper et al.l 
20041 : IProfumol 120051 : lAharonian et~aTI l2006aT) . Specifi- 
cally, we show that the improved angular resolution of 
CTA can differentiate these scenarios with surprising ac- 
curacy based on morphology alone, allowing for the con- 
struction of a rigorous model for the 7-ray emission from 
the GC center source. In turn, the resulting morphol- 
ogy observed by CTA will stand as a crucial ingredient 
in the understanding of high-energy emission from the 
entire GC region, including possibly the extraction of a 
dark matter annihilation signal. 

2. MODELS 

The morphology of the 7-ray emission from tt° decay is 
dominated by the distribution of Galactic gas. In order 
to produce a rigorous model which takes into account the 
full 3D mor phology of the target gas density, we employ 
the maps of lFerriera (|2012t ) which include not only contri- 
butions from a spherical diffuse halo, but also from struc- 
tures associated with the SNR Sgr A East, from belts of 
molecular clouds, and most importantly, from the high- 
density gaseous disk known as the circum- nuclear ring. 
This ring-shaped cloud of gas is located between 1-3 pc 
from the GC, inclined 20° with respect to the Galac- 
tic plane, and contains gas densities approximately two 
orders of magnitude larger than in the surrounding GC 
medium (jBecklin et al.l Il982h Bradford et al.l 120051) . In 
this work, we assume the central position and gas den- 
sity for each gas structure, and set the volume filling 
factor o f each structure to match those provided in Ta- 
ble 1 of iFerrierd (|2012t ). Given the high energy of the 
injected protons in 7-ray models, we ignore all informa- 
tion on the temperature distribution of Galactic g as, as 
it is inconsequential for 7-ray production. 

In their analysis of the combined GeV and TeV s pec- 
trum, recent work b y bothlChcrnvak ova et al.1 (|2011| ) and 
ILinden etP al. (2012) employed a diffusion constant tuned 
in order to provide a sharp transition between diffusive 
propagation in the GeV energy regime and rectilinear 
transport at TeV energies. This transition, which occurs 



at diffusion constants of approximately 1.4 x 10 29 cm 2 s _1 
for a diffusion zone of 10 pc, must be finely set in order to 
correctly explain the extremely soft emission spectrum at 
energies of approximately 100-500 GeV. Thus, a generic 
prediction of all scenarios which employ a single proton 
injection spectrum to correctly match both the GeV and 
TeV 7-ray emission spectrum, is cosmic ray propagation 
which is transitioning from the diffusive to the rectilin- 
ear regime at TeV energies. In the rectilinear regime, 
the cosmic-ray density falls as r~ 2 , while in the diffusive 
regime the cosmic-ray density falls instead as r _1 , due to 
the square dependence of the escape time on the size of 
the diffusion region. 

In either case, we will assume that cosmic ray pro- 
tons in this energy range tend to interact with gas less 
than one time before leaving the diffusion region, validat- 
ing the assumption that the energy spectrum of cosmic 
ray protons is pos ition independe nt. For the gas density 
maps provided by iFerrie rc ( 2012j) this assumption holds 
so long as the average value of the diffusion constant in 
the inner 10 pc of the galaxy exceeds 1.8 x 10 27 cm 2 s -1 . 
Under this restriction, coupled with the assumption that 
the 7-ray emission from the galactic center is in steady 
state, our models for the morphology of the 7-ray emis- 
sion depend only on the radial density of cosmic-ray pro- 
tons, rather than the exact diffusion scenario which pro- 
duces the proton morphology. 

Throughout this work, we will calculate the expected 
7-ray morphology for both the diffusive and rectilinear 
models of proton-propagation, finding that the qualita- 
tive arguments presented throughout this paper do not 
depend on the exact diffusion model considered. Fur- 
thermore, we can consider these two scenarios to bound 
the distributions expected for relativistic particle motion 
due to interactions with magneto-hydrodynamic waves, 
and can also be used as approximations of the expected 
effect in scenarios involving either non-homo genous dif- 
fusion scen arios such as those put forth by iFrver et all 
(|2007f) and lBallantvne et all (|2007l) . 

Finally, a significant uncertainty in this model pertains 
to the angular dependence of cosmic-ray injection from 
the GC point source: while a spherically sy mmetric dis- 
tribution was found by ILinden et al.l (|2012l ) to provide a 
compelling match to the current H.E.S.S. point source 
limit on the GC emission, other models are possible. We 
will comment on the assumption of isotropy in the cosmic 
ray distribution in Section [4] Throughout this work, we 
restrict our analysis to examining photons in the energy 
range of 1-10 TeV, which yields several simplifications to 
the analysis. First, the spectrum in this region is rela- 
tively flM J _Joll_owiu£abe^t fit spectral index T = 2.10 ± 
0.04 (|Aharonian et alJl2009| ). Secondly, the PSF of both 
H.E.S.S. and CTA are relatively constant in thi s region 
(|Aharonian et al.ll2006cl ICTA Consortiumll2011| ). 

The most up-to-date analysis of the GC with H.E.S.S. 
consisted of 93h of live-time with the instrument oper- 
ating in "Wobble" mode with a an average distance of 
0°.7 from the position of Sgr A *, and producing event s 
with a mean zenith ang l e of 23 ° (|Aharonian et al.1 l2009h 
From lAharonian et al.l (|2006bl Fig. 13) we infer an ef- 
fective area of 2 x 10 9 cm 2 , with only negligible varia- 
tion over the energy range 1-10 TeV. In this region, we 
adopt a point-spread function which is constant in energy 
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and foll ows the functional form given in lAharonian et al.l 
( 2006b) of a two-component Gaussian where the proba- 
bility of finding a photon in a radial bin d6 is given by: 



P(9) = A0(cxp(-^-2 ) + A c iexp(- 



(i) 
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Fig. 1. — Observed flux at the solar position as a function of the angle from the galactic center from the actual emission morphology (i.e. 
for a machine with a perfect angular resolution, top row), the emission as observed by an instrument with the angular resolution of the 
HESS telescope (middle row) and for an instrument with the angular resolution of the CTA telescope (bottom row) in the case of hadronic 
emission in the regime rectilinear proton propagation (left column), hadronic emission in the regime of diffusive proton propagation (middle 
column), and point source emission. The fluxes are shown logarithmically and binned to regions of 10 _4 °. The pointsource emission for 
an instrument with perfect angular resolution is not shown, as it would would provide a delta function at the center of the image. The flux 
shown stands as the flux between 1-10 TeV in energy averaged over 10000 realizations of 93 hour HESS observations and 1000 realizations 
of 500 hour CTA observations. 

2 x 10 10 cm 2 , which exceeds the H.E.S.S. effective area by 
an order of magnitude, and we adopt an equivalent func- 
tional form for the point-spread function as described in 
Eq. ([1]) for H.E.S.S., but rescale the parameters u\ and 
(72 such that the 68% containment radius of the photon 
signal is 0.03°. This yields cti = 0.0115 and a 2 = 0.03. 
While CTA contains additional improvements over cur- 
rent Cherenkov telescopes, especially stemming from its 
significantly lower energy threshold, the poorer angular 
resolution in the lower energy regime mitigates the effec- 
tiveness of CTA to test the morphology of the dense GC 
region. In this work we evaluate the performance of CTA 
after both 100 and 500 hours of observation, indicating 
both a conservative lower bound and a targeted obser- 
vation time for the GC region. This yields an expected 
14323 and 71613 photon counts, respectively. 

In order to simulate observations of the GC with both 
the H.E.S.S. and CTA instruments we employ Monte 
Carlo techniques to calculate the expected distribution 
of observed photons. We first calculate the 3D mor- 
phology for the true photon direction. In the case of 
hadronic emission, we calculate this by multiplying the 
3D distribution of gas with both the r -2 and r" 1 cosmic- 
ray densities assumed for rectilinear and diffusive trans- 
port of TeV protons. For dark matter models we assume 
a density distribution p{i) — r Q and evaluate scenar- 
ios a={1.0, 1.4}, as we note that current hydrodynam- 



with cti = 0.046, cr 2 = 0.12, A re i = 0.15 and A a nor- 
malization constant. However, in this work, these pa- 
rameters are set specifically to account for the spectral 
characteristics of the Crab Pulsar and a 60° zenith an- 
gle. These parameters provide a 68% c ontainment angle 
of 0.12 ° degrees. In the case of the GC, lAharonian et al.l 
(|2006d ) yields a 68% containment angle of 0.08°, and we 
thus linearly scale down the parameters 9\ and 02 to the 
values 0\ = 0.031 and 2 = 0.08 in order to obtain the cor- 
rect 68% c ontainment angle. Fin ally, H.E.S.S. observa- 
tions from ()Aharonian et al.1 ,2009) find a best fit intensity 
above 1 TeV oil >1Te v = (1-99 ± 0.09) x 10" 12 cm~ 2 s" 1 . 
Given the calculated effective area of H.E.S.S, this im- 
plies a point source observation of 1332 photons with 
energy in the range 1-10 TeV. 

In order to model the instrumental performance of 
CTA, we adopt best fitting pa rameters following the de - 
sign specifications set forth in lCTA Consortium! (|2011l ). 
noting, however, that the ultimate design specifications 
for the instrument are presently unknown. Specifically, 
we adopt an effective area in the 1-10 TeV band of 



4 



250 




3 
Z 



0.05 
Distance from GC ( 



- Point Source 

- Hadronic 



xVdof=1037 
V CTA 500h 




250 



4000 



-2000 



0.02 0.04 
Distance from GC (*) 

Fig. 2. — Expected photon counts as a function of the dis- 
tance from the GC (° from Sgr A*), for models where the 1-10 
TeV signal is generated by photons from a point source at the 
position of Sgr A* (black) or via hadronic emission from the cen- 
tral point source which propagates rectilinearly (p(r) oc r — 2 ) and 
subsequently interacts with gas (red). The mean shown in each 
model is the average of 1000 realizations, and error bars indicate 
the \J counts for a the average simulation. The value of x 2 /d.o.f 
is computed via a two-sample K-S test which does not depend on 
the binning used in the figure. 

ical simulations indicate the possibility that the inner 
dark matter de nsity profile is adiabatically contracted 
(i.e. a > 1.0) (iBlumenthal et all U985 iRvden fc Gmml 
[l987t iGnedin et al.l 12004 1201 llp These emission profiles 
are integrated over the line of sight, and then photons 
are selected from this distribution and smeared with the 
PSF of each instrument. We run 1000 simulations of all 
models in order to achieve reasonable statistical accu- 
racy, unless otherwise noted. 

3. RESULTS 

In Figure [T] we show 7-ray intensity maps depicting 
the actual emission morphology (i.e. for a machine with 
perfect angular resolution) (top row), the emission as 
observed by a machine with the angular resolution of 
the HESS telescope (middle row) and the emission as 
observed by a machine with the angular resolution of 
CTA (bottom row) in the case of hadronic emission in 
the regime rectilinear proton propagation (left column), 
hadronic emission in the regime of diffusive proton prop- 
agation (middle column), and point source emission. In 
models employing the effective area and angular resolu- 
tion of the H.E.S.S. telescope, the morphological features 
stemming from the gas density employed in the hadronic 
emission scenario occur on angular scales significantly 
smaller than the a = 0.08° angular resolution of the tele- 
scope blurring out the angular features which distinguish 
the hadronic scenario. Since these angular features are 
themselves centered around the position of Sgr A*, differ- 
entiating between the point-source and hadronic models 
becomes extremely difficult. 

Calculating the average cumulative-distribution func- 
tion over 1000 simulations of both a GC point source 
and the hadronic emission scenario, we employ a K-S test 
and find that current H.E.S.S. results would only be able 
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Fig. 3. — Same as Figure[2]but for protons which are propagating 
diffusively (p(r) oc r" 1 ) (red) 

to differentiate between models with a x 2 /d.o.f. = 0.85, 
which falls far short of providing a minimum 2tr level 
of confidence between the emission profiles. In Fig- 
ure [5] (top left) we provide a binned (at 0.01°) sample 
of photons collected by the H.E.S.S. telescope in each 
scenario. In realistic observations, this result is further 
complicated by residuals stemming from both the Galac- 
tic plane integrated over the line of sight, as well as from 
contaminating cosmic-ray backgrounds, both of which 
should appear isotropic in the small region under con- 
sideration. While additional H.E.S.S. observations time 
may slightly improve these statistics, these additional 
backgrounds make it unlikely that the H.E.S.S. telescope 
will be capable of differentiating between the point source 
and hadronic emission scenarios. 

In the case of CTA, the improved angular resolution 
will provide a much sharper view to distinguish between 
a point source and hadronic models. Furthermore, the 
greatly increased effective area will (in our simplified 
model where additional backgrounds are rejected) in- 
crease the x 2 mismatch between models linearly. Using 
100 hours of CTA observation, a K-S test provides a fit 
X 2 /dof = 208, which provides more than 14er differentia- 
tion between models. We find that the improved angular 
resolution of CTA will allow for a 3a rejection of the 
poorer fitting model with only ~5h of pointed observa- 
tion! In Figure [5] (top right), we again show a binned 
analysis for this dataset, noting specifically the under- 
density of photons observed within the inner 0.01°, which 
provides an independent, statistically significant indica- 
tion that would be difficult to explain with an additional 
diffuse or cosmic-ray background. We note that over the 
projected lifetime of CTA, nearly 500h of GC observation 
are expected, which would lead to a 32<r differentiation 
between models, with a result that is plotted in Figure [2] 
(bottom right). A wealth of information on the nature of 
the GC source will clearly be available in this case, going 
well beyond simply distinguishing between a point-source 
emission and a hadronic model. 

In Figure [3] we show the same observations as in Fig- 
ure [2] but under the assumption that high energy pro- 
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Fig. 4. — Cumulative Distribution Function (CDF) of photons 
observed within 0.05° of the GC for models of the total number of 
photons produced within ~0.07° of the GC. We show the case of a 
point source at the position of the GC (solid black), the hadronic 
model as described in Section [5] (green short dash), dark matter 
following an NFW profile a = 1.0 (blue long dash) and dark matter 
following a steeper profile with a = 1.4 (orange dot-dashed). 

tons travel diffusively through the inner galaxy, and 
thus follow a density distribution p(r) oc r -1 instead of 
p(r) oc r~ 2 . In this case, we find that the statistical differ- 
entiation between the Point Source and Hadronic scenar- 
ios slightly increases in all cases, due to the more diffuse 
nature of the energetic protons. However, the qualitative 
results are unchanged. Specifically, in the case of H.E.S.S 
observations this improvement is still insufficient to dif- 
ferentiate the two signals, as our K-S test obtains a fit 
X 2 /dof=1.49. 

In realistic models, additional emission sources must 
also be considered, including an isotropic cosmic-ray 
background, a line of sight background through the 
Galactic plan^3, and unresolved sources in the region sur- 
rounding the GC - all of which will contribute additional 
uncertainties to the differentiation of the point source 
and hadronic models. One particularly interesting back- 
ground could stem from the annihilation of dark mat- 
ter particles in the GC region. The morphology of the 
dark matter annihilation is partially constrained to be 
spherically symmetric with a density distribution which 
follows a form p(r) oc r~". While a standard value, 
a = 1.0 is employed in the st andard NFW dark mat- 
ter model (Na varro et al.|[l"997T) . the dark matter profile 
is highly uncertain in the GC region, and the gravita- 
tional effect from baryons in the GC may significantly 
steepen the dark matter dis t ribution ([Bl umcnthal et al.1 
[19M IRvden fc Gum] fl987t IGnedin et al.l I2004L 1201 i f. 
This makes it potentially difficult to differentiate be- 
tween dark matter models for TeV emission from the 
GC, and the possible combination of emission from both 
point source and hadronic sources. 

In Figure|4]we plot a projection for the cumulative dis- 
tribution function of photons observed with 500h at CTA 
as a function of the angular distance from the GC, for 

2 We note that this is approximately isotropic for the very small 
angular regions considered here. 
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Fig. 5. — Expected photon counts as a function of the azimuthal 
angle <j> = arctan(fe/Z) (shown in the inset), for both the point- 
source (black squares) and the hadronic emission scenario for pro- 
tons which propagate rectilinearly (red triangles) for photons ob- 
served within 0.05° of the GC. In the case of a point source, pho- 
tons are symmetrically distributed around the position of the GC, 
while the majority of emission due to hadronic injection follows the 
distribution of gas, which is aligned more closely with the Galactic 
plane. 

models of point source emission, hadronic emission, dark 
matter annihilation with a density profile a — 1.0 and 
with a density profile a — 1.4. We note that the point 
source model contains the fastest rising CDF possible, 
with a morphology uniquely determined by the instru- 
mental PSF. Any combination of emission from a point 
source and the hadronic model must produce a CDF 
which lies between the individual models, and the rel- 
ative contribution of each source class can be accurately 
(to within ~10% errors after 500h of observation) deter- 
mined by examining the CDF observed by CTA. How- 
ever, a small contribution from dark matter annihilation 
following an index a = 1.0 along with a dominant contri- 
bution from the point-source, may be misinterpreted as 
emission stemming from important contributions of both 
the point-source and hadronic models. Moreover, models 
where the emission is entirely dominated by dark matter 
which is highly peaked towards the GC (such as a — 1.4), 
may also be misinterpreted as some linear combination 
of hadronic and point source contributions. This uncer- 
tainty is a standard result from an attempt to identify 
three unknown intensities using only one constraint pa- 
rameter. 

A separate measurement is therefore necessary in or- 
der to constrain the relative contributions from all three 
source classes. An obvious choice is to model the angu- 
lar distribution of photons around the GC, noting that 
both the point source and dark matter models are spher- 
ically symmetric. In Figure [5] we plot the expected az- 
imuthal angular distribution for both the point source 
and hadronic models, counting the number of photons 
from a given angle <f> = arctan(6/^), i.e. the angle formed 
between the Galactic plane and the direction joining the 
GC and the photon location in the sky (see inset). We 
restrict the counts to photons within 0.05° of the GC, 
where contributions from the H.E.S.S. point source are 
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Fig. 6. — Same as Figure[5]but for protons which are propagating 
diffusively (p(r) oc r" 1 ) (red) 

believed to be largely dominant. While a (spherically- 
symmetric) point source provides a flat distribution in 
the incoming photon angle (as expected), contributions 
from the hadronic scenario deviate significantly, and are 
primarily aligned with the Galactic plane, due to contri- 
butions from both the circum-nuclear ring, as well as a 
contribution at angles of near 0° from hadronic interac- 
tions within the molecular cloud M-0.02-0.07. We note 
that with 500h of observation, an evaluation of the na- 
ture of the central TeV source can be made with more 
than ller confidence - without any reference to the radial 
distribution shown in Figure [21 Most importantly, this 
implies that in cases where the point-source, hadronic 
sources, and dark matter annihilation all contribute non- 
negligibly to the total TeV galactic center source, we 
can determine the relative contribution stemming from 
hadronic emission to within 20% with more than 2a con- 
fidence. In Figure [5] we show the same analysis as in 
Figure [51 under the assumption of diffusively propagat- 
ing protons and find the results to be qualitatively un- 
changed, with a x 2 /dof which improves by about 25%. 

4. DISCUSSION AND CONCLUSIONS 

We have shown that while current TeV instruments 
are incapable of distinguishing between point-source and 
hadronic emission models for the GC, the upcoming CTA 
will definitively differentiate between the signals within 
its first hours of observation at the GC, and will deter- 
mine the relative importance of each contributing source 
class at the 10% level over the course of its lifetime (see 
Fig. [Hand related discussion). These source classes can 
be distinguished based on either their radial or angular 
morphology, yielding two independent handles for the 
determination of the 7-ray source. The angular informa- 
tion is particularly important, as it allows the hadronic 
scenario to be easily separated from any combination of a 
GC point source and models of dark matter annihilation. 

While in this work, we have considered only 7-rays with 
an energy exceeding 1 TeV in our quantitative evaluation 
of the ability of CTA to distinguish between the hadronic 
and point-source signals, we note that the discrimina- 
tion power of CTA could be further extended depending 



on the specific instrumental point spread function which 
the instrument attains at lower energies. The continu- 
ation, or even softening of the E -20 7-ray spectrum to 
energies as low as 1 GeV has been demonstr ated by the 
combi nation of HESS (Aharonian et al. 2009) and Fermi- 
LAT (IVitale et all [200 91: IHooper fc Goodenougbl [20TTI: 
IHooper fc Linden! 120111 ) results.. This has the potential 
to create a much stronger statistical test, so long as the 
angular resolution of CTA remains on the same scale as 
the size of the circum-nuclear ring. Moreover, the differ- 
entiation between the Hadronic and point-source mod- 
els improves in the diffusive regime, which is thought to 
control the propagation of the protons which create the 
7-ray signal below ~100 GeV (jChernvakova et a"ill2011l ; 
ILinden et aT1l2012f) . 

We further note that while this study quantitatively 
examines only one specific scenario for hadronic diffu- 
sion through the GC region where the diffusion constant 
was fine tuned in order to match the softening of the 7- 
ray s pectrum at energies ~100 GeV (C hernvakova et al] 
120111 ). the model remains relevant for other physical 
models which have been proposed to control the propaga- 
tion of high energy protons in the galactic center region. 
Specifically, iFrver et al.l ([2007D employs a model which 
allows rectilinear propagation out to approximately 1 pc, 
after which the shocked winds produce a proton density 
which remains relatively flat between 1-3 pc, correspond- 
ing with the region dominated by the circum-nuclear 
ring. The comparison of our results in the case of both 
rectilinear (p(r) oc r~ 2 ) and diffusive (p(r) oc r _1 ) propa- 
gation show that the overall change in the pion produc- 
tion morphology is relatively unchanged by this factor of 
3 change in the proton density at the outer edge of the 
circum-nuclear ring. In fact, from our results, we would 
expect a marginal i mprovement in the differentiation be- 
tween the model of (jFrver et al.ll2007t ) and that of point- 
source 7-ra y production. The sam e is also true in the 
models of (Ballant vne et all I2007T) . where the diffusion 
constant varies inversely with the local gas density, pro- 
ducing a density of cosmic-rays which increases greatly 
in the region of the circum-nuclear ring. 

There is however, one important caveat to the scenario 
employed here. The morphology of observed high energy 
gamma-rays from the Hadronic scenario may become sig- 
nificantly less extended in cases where high energy pro- 
tons are not isotropized by the galactic medium through 
which they are propagating. In this case, 7-rays observed 
at the solar position will preferentially stem from protons 
which were themselves originally pointed towards the so- 
lar position, due to the relativistic beaming of pions pro- 
duced in p-p scattering. This complication may become 
important in scenarios where we have assumed rectilinear 
proton propagation, which in the case of a simple, homo- 
geneous diffusion constant, would imply that protons do 
not become completely isotropized before leaving the dif- 
fusion zone. In this scenario, we find that for diffusion 
constants exceeding 9.3 x 10 29 cm 2 s -1 , the random walk 
approximation for particle diffusion predicts that the av- 
erage proton will not undergo a change in direction be- 
fore leaving the 10 pc dif fusion region. In the stud ies of 
iChernvakova et all ([20111 ) and ILinden et all (|2012[ ) this 
corresponds to a proton energy of only 2.2 TeV. It is these 
protons which dominate the 1 TeV signal, since the E~ 2 
proton injection signal implies that the most important 
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contributor to 7-rays of a given energy are protons of a 
very similar energy. This could mean that some fraction 
of protons which are undergoing p-p scattering in our 
simulations are not isotropically distributed. 

However, the random-walk approximation for proton 
diffusion is notably poor in this region, and all protons 
will undergo some scattering off of the magnetic field 
inhomogenities in the diffusion region. The exact de- 
gree of proton anisotropy will depend sensitively on the 
specifics of this diffusion scenario, and specifically on any 
inhomogenities in the magnetic field structure, which are 
both invariably present, and difficult to directly deter- 
mine. Lastly, we note that even if the vast majority of the 
signal is anisotropic (thus appearing very similar to the 
point source calculation), even a small isotropic compo- 
nent could be distinguished from the point source calcu- 
lation, given the incredible x 2 detection of the Hadronic 
model by the CTA. 

We note that this same process could also place con- 
straints on the dark matter annihilation cross-section 
- using the observed angular and radial profiles to de- 
termine the maximum contribution of photons following 
a morphology consistent with dark matter annihilation. 
However, due to the small angular region considered, as 
well as the extremely bright point-source emission ob- 
served by H.E.S.S., these constraints are not competitive 
with those determined, for example, by H.E.S.S. observa- 
tions o f the regions directly abo ve and below the Galactic 
plane (jAbramowski et alJl201lh . While the GC remains 
an extremely interesting region for setting dark matter 
constraints with CTA, the best limits will likely continue 



to be set by analysis of regions directly off of the Galactic 
plane. This scenario may become interesting, however, 
in cases where the dark matter density in the GC re- 
gion is found to be highly adiabatically contracted (e.g. 
a > 1.7), leading to extremely enhanced fluxes within 
the inner pc of the Galaxy. 

Pinpointing the nature of the GC source will allow 
us to more effectively search for a dark matter signal 
in this region, as the differentiation between the point 
source and hadronic models would allow for an extrapo- 
lation of their expected emission profiles into the regions 
where strong dark matter limits can be set. Addition- 
ally, the extrapolation of the observed CTA emission 
will greatly refine Fermi-LAT models of the morphol- 
ogy of the central point source, as shown in lLinden et al.l 
(2012). Thus TeV 7-ray observations will be critical to 
validate or constrain particle dark matter models that 
could explain other observations tentatively indicating 
sign als from dark matter annihilation at lower energies 
(e.g iHoopen (|2012f) ). In this scenario, the role of CTA 
will be highly complementary to an extended Fermi-LAT 
campaign, especially in accurately subtracting a diffuse 
signal emission from the bright, and possibly extended, 
central source. 
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